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Abstract. We propose that the first stage of membrane fusion need not be the formation of a single 
stalk. Instead, we consider a scenario for defect-induced membrane fusion that proceeds cooperatively via 
multiple stalk formation. The defects (stalks or pores) attract each other via membrane-mediated capillary 
interactions that result in a condensation transition of the defects. The resulting dense phase of stalks 
corresponds to the so-called fusion intermediate. 
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When two bilayer membranes approach each other, 
they may fuse to form a single bilayer membrane. This 
process (and its reverse) are of great importance for many 
processes in a living cell. Nevertheless, the mechanism by 
which membranes fuse is still a matter of debate PI21I31- 

The most widely used description of membrane fu- 
sion assumes that it involves several steps. It has been 
argued 4.5 6 7.8 9. 1 0II11II12| that the initial connection 
between the membranes is formed either through a stalk 
or through a pore |7lll2j . In fact, it is now commonly ac- 
cepted that the initial inter-membrane contact is, most, 
likely, a stalk [2j (see cartoon Fig.^ a )- This initial state 
is called the hemifusion A recent theoretical anal- 
ysis of the free-energy cost of the hemifusion state [TT)| 
predicts that, depending on the magnitude of the spon- 
taneous splay of the lipids, this free energy F can either 
■ be positive (F ~ 45 fcs T for the common case of DOPC 
; lipids) or negative (F ps -30k B T for DOPE lipids that 
have a large and negative spontaneous splay). 

The second step of the fusion reaction is the expan- 
sion of the hemifusion zone (diaphragm) 1 1 . The analysis 
of ref. 55 predicts that the expansion of the hemifusion 
diaphragm is energetically favorable only if the sponta- 
neous splay modulus of lipids is large and negative (e.g., 
for DOPE lipids) while it costs energy for membranes com- 
posed of other types of lipids (e.g., DOPC). It is believed 
that in the case when the expansion of the hemifusion di- 
aphragm is energetically unfavorable, fusion proteins gen- 
erate the additional driving force needed to expand the 
diaphragm. The fusion process is then completed by the 
subsequent formation and expansion of a fusion pore. 

A direct experimental conformation of the above sce- 
nario is still lacking after more than two decades of in- 
vestigations. Only recently [141113] Yang and Huang have 
succeeded to crystallize a stable phase of membrane stalks 
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in multi-lamellar system, and verify the predicted struc- 
ture of the stalk intermediate. No similar experiments have 
been reported in the case of two membranes. 

Our principal hypothesis in this paper is based on the 
observations that (i) two stalks or two pores attract each 
other, (ii) the translational entropy associated with the 
formation of a single stalk or pore is sufficiently large to 
allow the spontaneous emergence of a dilute gas of such 
defects, even if the elastic free energy of a defect is positive. 

We argue below that the defects that are thus formed, 
will attract each other and self assemble into a struc- 
ture that has the characteristics of a hemifusion zone (see 
Fig. He). This is a cooperative effect, similar to conven- 
tional order-disorder transition. We emphasize that the 
predicted aggregated phase of stalks constitutes an inter- 
mediate and not a final stage of membrane fusion. The 
final stage can proceed e.g., either through stalk coales- 
cence and a subsequent pore formation, or through the 
expulsion of stalks via budding. This final stage of fusion 
is beyond the subject of the present work. 

We stress that the assumption of multiple stalk for- 
mation is not as farfetched as it may seem. In fact, the 
transitions from the lamellar to the inverted hexagonal 
phase |1KU17| . and from the lamellar to sponge phase [T%] 
are both examples of a similar effect, where the lamellar 
phase is transformed into a highly connected structure. 
In any event, the proposed mechanism for the formation 
of the intermediate state in the fusion process is not re- 
stricted to the case of spontaneous (passive) multiple stalk 
formation. Our conclusions also apply to the case where 
fusion proteins actively facilitate stalk formation [TTI19I 

EDI- 

To model the condensation of defects during mem- 
brane fusion, we exploit the close analogy between the 
current process and the interaction between mobile re- 
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ceptors that arc responsible for the adhesion between cell 

In what follows, we consider the scenario proposed by Bru- 
insma, Goulian, and Pincus in ref. |22| . 






Fig. 1. Schematic representation of the cooperative forma- 
tion of the fusion intermediate state via condensation of the 
self-assembled defects (stalks). The stalks represent elemen- 
tary defects which interact via an effective, membrane-induced 
mechanism, (a, 6): Two bilayer membranes with two interact- 
ing stalks. (6): Within the framework of the model, the mem- 
branes have an equilibrium separation H. The stalks constrain 
undulations of the membranes. The membranes overshoot the 
equilibrium separation H in the vicinity of the stalks. This is 
the driving factor for the attraction between the stalks. When 
the stalks are close to each other, they can share the com- 
mon overshoot and thus reduce the total bending energy, (c): 
Cartoon of the two-dimensional (2D), phase-separated (con- 
densed) phase of stalks in the contact zone of two membranes. 
The dense phase of stalks is equivalent to a 2D liquid of defects. 
The stalks are not the only type of possible defects - pores, or 
pores coexisting with stalks can be also formed. 



Consider two bilayer membranes separated by an equi- 
librium distance H. This equilibrium separation can be 
a minimum of the inter-membrane interaction potential 
V(h) |35| . where h is the inter-membrane separation, or 
an optimal distance maintained by an external force. In 
the case of biological fusion, H is determined by inter- 
membrane proteins and by the glycocalyx coating (see ref. 

for a review). We are interested in the self-assembly 



of junctions (i.e., stalks or/and pores) between the mem- 
branes (see Fig. Each junction represents an elemen- 
tary defect. 

We assume that each defect imposes a local membrane 
spacing Hq that is different from H . In the present case 
of fusion, Hq is simply zero, as implied by the geometry 
of a stalk, or a pore |l()llllj . This imposes the boundary 
conditions on the inter-membrane spacing at the point of 
a junction po [22]: h(po) = 0, and V/i(po) = 0. The free 
energy of the system can then be written in the form |22| : 



F 



d 2 p 



f(V 2 M 2 



V"(H) 



{h-Hf 



(1) 



where the first term is the Helfrich bending energy 
with the bending modulus k, the second term is the inter- 
action energy between the membranes; h(p) depends on 

the lateral coordinates p = (x,y), and V"(H) = — gf ^ ' 
[the deviations of the inter-membrane separation h(p) from 
the equilibrium value H are assumed to be small, and 
thus the inter-membrane interaction potential V(h) is ex- 
panded to quadratic order in (h — H)]. 
The minimum of Eq. Q is given by: 



V 4 h 



h-H 







(2) 



where A = [k/V" (H)] 1 / 4 is the capillary length - the char- 
acteristic length of the perturbation decay. This length is 
analogous to the stalk width R in the recent model of the 
stalk by Kozlovsky and Kozlov |TJ]]. The membrane profile 
modified by the presence of a junction located at p = 
is given by the solution of Eq. J2Jl with the corresponding 
boundary conditions, h(0) = 0, and Vfr(0) = [22]: 



h(p) =H+-Hkei(p/X), 

7T 



(3) 



where kei(a;) is the Kelvin function. An interesting prop- 
erty of this profile is that it overshoots the equilibrium 
inter-membrane separation, H. This effect is even more 
pronounced if the non-linear contribution to the interac- 
tion potential is taken into account |37II39| . The effect of 
strong overshooting of membrane profiles when pinched 
together by optical tweezers has been observed experimen- 
tally by Bar-Ziv et al. [23 and analyzed theoretically by 
Menes, Safran, and Kessler 



The free energy of a single junction is obtained by sub- 
stituting the membrane profile h(p) into the free energy 
Eq. <[1]>: 



F = 4k- 



IP 



(4) 



This free energy can be directly related to the free energy 
of the stalk ^H]- The dimensionless parameter \ = tt can 
thus be obtained for a given set of the inter-membrane dis- 
tance H and the stalk width A (R in the notations of ref. 
|1(J|). For example, the free energy of the unconstrained 
stalk in the case of DOPC lipids was estimated in ref. ^U] 
to be F ss 43/cbT. This gives x ~ 1> f° r a typical value 
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of the bending modulus k sa lOfc^T. This implies that 
H/X = 0(1). This finding is consistent with the values for 
H and A which follow directly from the calculation of the 
stalk profile in ref. jT|]]: H « 6.2 nm and A « 8.7 nm, lead- 
ing to the estimate \ ~ 0-7- In other words, the present, 
simple model allows us to interpret the results of the so- 
phisticated computation of the stalk energy of ref. |1()| in 
terms of a single parameter that follows directly from 
the two intrinsic parameters (H and A) of the stalk model. 
In fact, it is argued in ref. ^Oj that the stalk energy is 
minimal when the stalk width A is of the order of the 
inter-membrane separation H. With that information, it 
follows that the free energy of a single stalk F w 4 k. 

We stress that within our model, a stalk can only 
have a positive, elastic energy. This is because we assume 
that the only effect of a stalk is the constraint on the 
inter-membrane separation, and we neglect the topologi- 
cal change upon the stalk formation (see Fig.Q. To take 
this into account, each monolayer of membranes must be 
treated separately 4(5]. This topological change can lead 
to a negative stalk energy for lipids with sufficiently large 
and negative spontaneous splay ^U]. Within the frame- 
work of our model this can be taken into account phe- 
nomenologically by considering the stalk free energy e s 
as a sum of the elastic contribution F and the core free 
energy F core : e s = F + F core , where F core can be either 
positive or negative. 

The existence of the attraction between the junctions 
follows from the analysis of the free energy of a lattice 
of junctions |22- In the limit p ^> A, the interaction free 
energy between two junctions has the form |41|: 

where, C = « 25.5. The larger the bending rigid- 

ity k, the stronger the effective attraction, and the longer 
the range of this attraction. The characteristic range of 
the interaction, A is of the order of 10-20 nm for an ar- 
tificial phospholipid membrane under physiological con- 
ditions 23 29] and thus the long-distance limit Eq. 10 
should be accurate if the junctions are separated by more 
than 20 nm. We stress that in the opposite limit p <C A, 
the effective interaction remains attractive [52]. Hence, all 
qualitative conclusions hold irrespective of the value of p. 

To estimate the phase behavior of the stalks, we use 
a simple mean-field model analogous to the ones used to 
describe the aggregation of adhesion molecules or patches 
within the adhesion zone of two biological or biomimetic 

The mean-field free energy of the self-assembling junc- 
tions can be constructed with the effective inter-junction 
potential, Eq. JSJ, contributing to the second virial coeffi- 
cient: 

f = k B T(j) In + (1-0) In (1-0) (6) 
+ J0(l-0)+e s 0, 

where / is the free energy per elementary surface cell, 
(p is the surface fraction of the junctions, e s is the stalk 



free energy (e s acts as a chemical potential of defects), 
and J is the effective, thermodynamic interaction poten- 
tial between the junctions. J is obtained on the level of 
the linearized second virial coefficient: 

J = Cxk^, (7) 

where d = 16 \/2 tt sin(37r/8) « 65.7. Note that the 
coefficient C\ is remarkably large. Taking into account 
that a typical value of the membrane bending modulus is 
k ~ 10 — 20 ks T, the large value of C\ implies that the on- 
set of the phase separation of defects occurs at a very low 
value of the fusion control parameter, \ = -p-. Indeed, it 
follows from the analysis of / that the critical point of the 
"liquid-gas" phase separation of junctions is J c /ksT = 2 
and (f> c = 0.5. This implies, if we adopt k ~ lO/c^T, that 
the onset of the phase separation occurs at \c ~ 0.003. 
Comparing \c with the value of \ ~ 1> estimated above, 
one concludes that the phase separation occurs already 
at a vanishingly small concentration of defects. The im- 
portant message is that the higher the stalk energy, the 
higher the strength of the effective, inter-junction attrac- 
tion, and thus the smaller concentration of defects induces 
the phase separation. 

The above arguments indicate that stalk condensation 
should be possible for reasonable values of the param- 
eters characterizing biological membranes. Yet, the key 
question is: does it happen in practice? In fact, there 
is experimental evidence that supports the present sce- 
nario: very recently, Yang and Huang |14U15j reported 
x-ray scattering experiments that show the spontaneous 
formation of an ordered, dense multiple-stalk structure 
in a multi-lamellar system of bilayer, diphytanoyl phos- 
phatidylcholine (DPhPC) lipid membranes. From their 
scattering data, Yang and Huang were able to reconstruct 
both the global multiple-stalk structure of the fusing mem- 
branes, and the structure of the individual stalks. The 
latter shape turned out to be consistent with the classi- 
cal stalk structure [Tr)llTTir7lFT5] . We suggest that this kind 
of multiple-stalk structure should be present in practically 
all fusion experiments with artificial membranes, provided 
that the membranes or vesicles have large enough area of 
inter-bilayer contact (see e.g., ref. 

There are also recent computer simulations studies 031 
144) . that report either a simultaneous formation of two 
adjacent fusion zones [H] even in a fairly small fusing 
vesicles, or a coexistence of fusion stalks and pores 43 . 

In summary, we propose a possible mechanism of mem- 
brane fusion, via a multiple stalk formation. We predict 
that the intermediate stage of membrane fusion represents 
a phase-separated phase of self-assembled defects (stalks 
or pores, or both stalks and pores). Multiple defect gen- 
eration is a mechanism alternative to a hypothesis of a 
single stalk or pore formation and subsequent expansion 
of the fusion diaphragm. The physical origin of the pro- 
posed mechanism is the membrane-induced, effective at- 
traction between fusion defects, similar to the attraction 
between adhesion receptors. Our conclusions apply both 



4 



D. B. Lukatsky and Daan Frenkel: Multiple stalk formation as a pathway of defect-induced membrane fusion 



in the case of spontaneous stalk formation and in the case 
of protein-assisted stalk formation. 
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